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We show that the calculation of L-loop Feynman integrals 
in D dimensions can be reduced to a series of matrix multi- 
plications in D X L dimensions. This leads to a new type of 
expansions for the critical exponents in three dimensions in 
which all coefficients are calculable exactly. 

1. Recently, an efficient method has been developed for 
constructing all Feynman diagrams of a quantum field 
theory recursively to arbitrary order in perturbation the- 
ory [0. The construction has also been automatized by 
computer algebra [Q, so that one can now easily list all 
diagrams to any desired order, limited only by computer 
time. 

The only open problem is an efficient evaluation of the 
associated Feynman integrals. The highest loop number 
so far has been treated in ^''-theory where renormaliza- 
tion constants of field normalization and mass are now 
known numerically up to seven loops 

For a reliable calculation of all critical exponents with 
an accuracy which has become recently available in satel- 
lite experiments on superfluid helium Q it is desirable to 
know the perturbation expansions to even higher order 
1^,^. In this note we would like to propose a pertur- 
bative evaluation of many- loop Feynman integrals which 
reduces the integrations over D x L loop momenta to a 
series of multiplications of x L-dimensional matrices, 
thus reducing greatly the number of numerical opera- 
tions. 

2. The method is most easily introduced by treating an 
exactly solvable class of integrals 



D 



1 



(27r)^ (H-p2)«^ 



whose integration via Schwinger's proper-time technique 
yields 



1 Tia-D/2) 
(47r)^/2 f(a) ■ 



(2) 



Our perturbative evaluation starts out by introducing 
into (P a smallness parameter e, and considering the 
functions of e 
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{2Tr)° (1 + ep2)'^A 



(3) 



For this function, we can easily find a perturbation ex- 
pansion in powers of e, by rewriting it as follows: 



1 



d^p 



(47ra)^/2 J {TT/a)D/ 



,-Qp2-S.(p) 



with 



5,(p) = -log(l + ep^)-ap^ 



|(p^)^-y(p^)^ + j(pY 



(4) 

(5) 
(6) 



Considering ap^ as a free "action" , and S'e (p) as an "in- 
teraction" , we introduce the Gaussian expectation values 



(/(P)) 



re^^'P /(P), 



{Tr/a)D/2 

and rewrite the integral (jj) as a harmonic expectation 

1 



(7) 



(47ra)^/2 



,-5'.(p)\ 



(8) 



The exponential e"'^^^^-* is now expanded in a Taylor se- 
ries, and the expectation values {Si{p) ■ ■ ■ 5'e(p)) are cal- 
culated via a Wick expansion using the basic correlation 
function 

(P^iP^a) = Gji^a = 7^(5j,j2, i = l,...,D. (9) 
2a 



(1) The nonzero expectation values of higher products 



are calculated from the recursion relation 



Cjj...i2„ — ^ili2^i3U---i2re-li2n "I" ^3 ^^2 U ■ ■ ■i2n 

+ . . . + Giji2„-l^j2*3---j2n "I" ^ili2n^i2i3---j2n-l • 



(10) 



(11) 



For example, we find 
1 



(47r)l/2 
1 



3e 25e2 lOSe^ 16596"* 
^ ¥ ^ 128 ^ 1024 32768 ^ ' 

15e 385e2 4725e3 228459e4 



(87r)3/2 y - ■ 16 ' 512 ■ 8192 ' 524288 

The corresponding cumulant expansions are 



1 



(47r)i/2 
1 

(87r)3/2 



g 16 ^ 16 ^512*^ ^512*^ 



(12) 
(13) 



Evaluating the exponents at e = 1 successively as par- 
tial sums of orders 1,2,3,. . . ,15 [not all terms are written 
down in (|l^) and (|l^) , for lack of space] , the relative de- 
viations from the respective exact values 1/2 and I/Stt 
[which follow from Eq. (^)] are shown in Table | [columns 
marked with (ps)]. The deviation sequences tend slowly 
toward zero, and grow again for larger order as a conse- 
quence of the asymptotic nature of the expansions. 

The partial sums yield integrals whose first ten terms 
tend exponentially fast toward the exact results, as shown 
in Figs. |l|. 
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FIG. 1. Exponentially fast convergence of the first ten ap- 
proximations, obtained from the successive partial sums of 
the exponents in ( |l2[ ) and (^), against their respective exact 
values 1/2 and l/8n ~ 0.397884 ... . The number n denotes 
the order of the approximations. 

The convergence is accelerated by approximating the 
successive exponents a la Fade, leading rapidly to accu- 
rate results, as we can see in the right-hand columns of 
Table |. 

3. As a nontrivial example consider now the Feynman 
integral 



{2n 
1 



{2tt)D J {2n)D 
1 1 



Xe~[Pl+P2+P3 + (Pl+P2+P.3)^]-'S5(pi,p2,P3) ^-^^g-j 

with a free "action" Pi + P2 + P3 + (Pi + P2 + Ps)^, and 
an "interaction" 



'S'e(pi,P2,P3) = 



log(l + ep2) ^ p2 



log 



^Pa 

\a=l / 



(16) 



The free action is now rewritten as a quadratic form 
^PaiKai,bjPbj, with the matrix 



ai.bj 



2 1 1 

fc,6 = 2 ( 1 2 1 I , a = 1,2,3. 
1 1 2 



where 6ij is the Z3-dimensional Kronecker symbol. We 
now define harmonic expectation values 



(/(pi,P2,P3)>^(27r)3^/2(detif)i/2 
d^Pi fdDp2 fd^p3 



{2n)D J {2tt)D J {2tt)D 
where 

det K = det = 32^ 
The basic correlation functions are 



gab — fcab ^ o 



(17) 

/(pi,P2,P3)e-^P"^"-^^'"'^- 



ai,bj 


— Qab^ij ; 
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In terms of (^_7|) , the integral ( |15| ) can now be considered 
as a harmonic expectation value 

/^(e) = (27r)-^^/2(jg^^)-i/2^g-S.(pi,P2,P3)^^ (18) 

which is evaluated perturbatively by expanding the ex- 
ponential into a power series, and forming all Wick con- 
tractions. The results are for D = 1 and 2: 



1 



(27r)3/2 V32 
1 1 



1 27 I 9 2 , 24 3 I 249 4 
g32''^64'^ '4096 ^ 8192 ^ 



(27r)3 V32^ 



(19) 



l + p2 l+p2 l+p2 l + (pi+P2+P3)2 



(14) 



We introduce again the smallness parameter e as in (|^), 
and rewrite the integral (^4[) as 



jD,. fd^ fd^ fd^ 
^ ' J (2^)^ J {27r)D J {2n) 



D 



and the successive evaluation of the exponent at e = 1 
yields the approximations 

0.026097, 0.030038, 0.030236, 0.031169, (20) 
0.001195, 0.002304, 0.002823, 0.003356, (21) 

which converge against the precise values of the Feynman 
integral 1/32 = 0.03125 and 0.00424027 . . . , respectively. 
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These tend exponentially fast against the exact value, as 
seen in Fig. |[ 
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FIG. 2. Exponentially fast convergence of the approxima- 
tion sequences (^o|) and (^l|) against their respective exact 
values 1/32 and 0.00424027 .... The number n denotes the 
order of the approximations. 

The same Fade approximants as before yield the num- 
bers 



0.030894, 
0.003019, 



0.030171, 
0.003091, 



0.03064, (22) 
0.00309, (23) 



which, surprisingly, show no improvement with respect 
to the sequences (|20|), (|2l|). 

4. So far we have calculated only vacuum diagrams. 
There is no problem in extending the method to dia- 
grams with external legs. The first nontrivial Feynman 
integral in the self-energy of a 0^-theory, for example, is 



' ^^^^7 (^7 



(27r)^ 
1 1 



1 



i + (pi + q)2 i + pi i + (pi + P2)2' 

which we rewrite as 



(24) 



/^(q;e) = e-"^ , , / . xn 

' ' {2tt)0 J (2Tr)^ 

Xe^[Pl+P2 + (Pl+P2)^]-'5,(pi,P2,q)-2piq^ ^25) 

with a free "action" Pi + p^ + (pi + P2)^, a q-dependent 
"interaction" 



'5'e(pi,P2,q) = 



- log(l + epl) ~ pI 
e 



-log 
e 



I] Pa 



I] Pa 



(26) 



where we have set p3 = pi -I- p2. After removing the 
extra source term 2piq by a shift of the integration vari- 
able Pi ^ Pi — q, which produces an exponential factor 
g2feii(3 ^ integral can be treated as before, yielding an 
expansion in powers of e with q^-dependent coefficients. 

5. The above treatment of Feynman diagrams offers the 
interesting possibility of deriving a completely new type 
of power series expansions for the critical exponents of 
0'*-theories in three dimensions, in which expansion co- 
efficients are calculable exactly, in contrast to the known 
expansions of Ref . [|| . For this we observe that our power 
series for the Feynman integrals can be reexpanded in 
powers of g by inserting e = g/ (m + g). Then ll{e) and 
/Ke) become 



Ilia) 



^expfM_5l. 
(4^)1/2 c'^P 1^ 8 4 



1 



(47r)3/2 



exp 



15s 
16 ' 



lir 9 
64 8 
2355^ 
512 



23^3 
64 



where g = g/m. The coefficients have now alternating 
signs. These expansions must be evaluated in the strong- 
coupling limit g ^ oo. This is again possible with the 
help of diagonal Fade approximations, and the results 
are precisely the same as those in Table |. Note that the 
choice of the parameter m is irrelevant in this limit. 

Now we recall that the critical behavior of (/)^-theory 
can also be derived from power series expansions of the 
renormalization constants in the strong-coupling limit 
gs/niB oo, where gs and tob are the bare coupling 
strength and mass of the theory. Let us identify the above 
expansion parameter g with cgs/mB, where c is an ar- 
bitrary parameter which may be chosen at the end to 
optimize the speed of convergence. Then we insert the 
expansions of the Feynman integrals in powers oi gs/mB 
into the Feynman diagrams of the perturbation expan- 
sions of the renormalization constants of the 0^-theory. 
This leads to a new type of series in powers of gn/'m-B 
whose coefficients are determined exactly. By going to 
the strong-coupling limit gs/n^B — >■ oo as described in 
Rcf. 1^,^, we can obtain the critical exponents in three 
dimensions in a completely new way. 

From the connected vacuum diagrams we obtain an 
expansion of the form 



1 



aigB 



0255 



+ ...) 



(27) 



For a comparison with experiment one assumes the bare 
mass to go linearly to zero as the temperature T ap- 
proaches the critical temperature, i.e. cx (T/Tc — 1). 
The vacuum energy W has then for small m b the scaling 
behavior {T/T^-lf'"' oc m^^^a, where a = 2-3t/ is the 
critical exponent of the peak in the specific heat in three 
dimensions {via the exponent of the correlation length). 
Since gB grows for small niB like I/jtib, the power series 
for Z in ( p7| ) must behave like Z oc g^~^ ■ The critical 



3 



exponent a is therefore found from the strong-couphng 
limit of the power series 



dlogZ/dloggs), 



(28) 



evaluated as explained in Ref. 

Also the exponent co governing the approach to scal- 
ing can be obtained from the series ([28[), since ai^gs) 
approaches a for large gs like a+const/g^. Thus lo 
is found from the strong-coupling limit of the series 
^{gs) = -1 - d\oga'{gB)/d\oggBa'{gB) (see |§ for 
more details). 

The third independent critical exponent may be de- 
rived from the perturbative corrections to the field nor- 
malization. For this we take the expansion for the self- 
energy S(q) and extract the coefficient of the term. 
This is done most simply by differentiating all Feynman 
integrals of the two-point function using the rule (in D 
dimensions) 



d 



1 



92 



2D dq^dqfj 



and setting q = 0. The momentum q always flows 
through one or more lines of the diagram, allowing the 
differentiations to be performed as follows: 



52 



1 



-2D 



dq^dq,, (p - q)2 + m'^ [(p - q) 



,212 



+ 



i(p- 



[(P- 



(30) 



After setting q = in the resulting Feynman integrals, 
the calculation of the field normalization proceeds using 
the above formulas for the vacuum diagrams. 

Note that the e-dependent Feynman integrals are all 
finite for small e. However, as e increasing towards unity, 
a pole is encountered at e = 2/3. Thus, the expansion of 
a pole term cx l/(e— 2/3) should be subtracted from the e- 
expansions of divergent integrals before one can evaluate 
perturbation series at e = 1, i.e., before going to the 
strong-coupling limit of . 
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TABLE L Relative deviations of approximations 11(1) and 
11(1) from their exact values, once from successive partial 



(29) sums (ps), once from successive Pade approximations. 



n 


7^(1) ps 


/|(1) ps 


[N,M] 


Pade 


/|(1) Pade 


1 


0.17910 


0.49063 


[0,1] 


-0.391430 


-3.85230 




2 


0.06980 


0.30377 


[1,1] 


0.009816 


0.18603 




3 


0.02516 


0.19394 


[1,2] 


0.002101 


0.09851 




4 


0.00981 


0.12674 


[2,2] 


0.000935 


0.02490 




5 


0.00516 


0.08416 


[3,3] 


0.001821 


0.00309 




6 


0.00256 


0.05655 


[4,4] 


0.000637 


-0.00062 




7 


0.00027 


0.03836 


[5,5] 


0.000652 


-0.00062 




8 


-0.00021 


0.02625 


[6,6] 


0.000325 


4.8349 10' 


-5 


9 


0.00124 


0.01810 


[7,7] 


0.000325 


8.0048 10" 


-6 


10 


0.00114 


0.01255 


[8,8] 


0.000192 


3.3873 10" 


-5 


11 


-0.00273 


0.00875 


[9,9] 


0.000194 


1.6181 10" 


-5 


12 


-0.00275 


0.00613 


[10,10] 


0.000127 


1.8694 10" 


-5 


13 


0.01001 


0.00432 


[11,11] 


0.000130 


1.0960 10" 


-5 


14 


0.01000 


0.00305 


[12,12] 


0.000090 


1.1146 10" 


-5 


15 


-0.05026 


0.00216 


[13,13] 


0.000095 


7.3635 10" 


-6 



6. Summarizing we have shown that just as the genera- 
tion of all Feynman diagrams |Q] , also their evaluation can 
now be automatized by computer algebra, with the ad- 
ditional advantage that all expansion coefficients emerge 
as exact numbers. 

The actual calculation of diagrams with eight and and 
more loops will, of course, require a considerable amount 
of work. It remains to be seen whether the new perturba- 
tion expansions converge more rapidly against the exact 
results than the presently known purely numerical ones. 
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